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We show that respecting the underlying Lorentz structure in the parton model has very strong 
consequences. Failure to insist on the correct Lorentz covariance is responsible for the existence of 
contradictory results in the literature for the polarized structure function g2(x), whereas with the 
. . . correct imposition we are able to derive the Wandzura-Wilczek relation for g2{x) and the target- 

mass corrections for polarized deep inelastic scattering without recourse to the operator product 
. expansion. We comment briefly on the problem of threshold behaviour in the presence of target- 

' mass corrections. Careful attention to the Lorentz structure has also profound implications for 

the structure of the transverse momentum dependent parton densities often used in parton model 
treatments of hadron production, allowing the fcr dependence to be derived explicitly. It also leads 
^ ' to stronger positivity and Soffer-type bounds than usually utilized for the coUinear densities. 

o 

m 

^ ■ In this paper we shall show that imposition of the correct Lorentz structure in the parton model allows us to relate 
, two issues which, a priori, do not seem to be connected with each other: the derivation of a consistent expression for 
'""'i the polarized deep inelastic scattering (DIS) structure function .92(2^) in the parton model, and the derivation of the 
I \ higher twist target-mass corrections i.e. corrections of the form M'^ /Q"^, where M is the nucleon mass, to gi{x,Q^) 
j> . and g2{x,Q'^). 

' As will be discussed below, the target-mass corrections have previously been derived in a very complicated way from 
I/"} the operator product expansion (OPE), and the derivation of the Wandzura-Wilczek (WW) expression for g2{x) has 
involved a questionable analytic continuation in the OPE moments jj. It is thus particularly interesting that these 
, results can be derived in a field theoretic context without use of the OPE. Target-mass corrections for unpolarized 
^\ • DIS were first derived by Nachtmann employing a very elegant mathematical approach in which the power series 
' expansion used in the OPE was replaced by an expansion into a series of hyperspherical functions (representation 
0^ functions of the homogeneous Lorentz group). Later, also within the context of the OPE, Georgi and Politzer Q 
re-derived Nachtmann's results using what they called an alternative analysis "for simple-minded souls like ourselves" 
i.e. based on a straightforward power series expansion but, in fact, requiring a very clever handling of the combinatoric 
aspects of the problem. 

^ , The derivation of target-mass corrections for polarized DIS turned out to be much more difficult. Several papers 
' 0i Q succeeded in expressing the reduced matrix elements a„, dn of the relevant operators in terms of combinations 
of moments of the structure functions, but did not manage to derive closed expressions for the structure functions 
gi^2 themselves. The latter was finally achieved in 1997 by Piccione and Ridolfi ^] and later generalized to weak 
interaction, charged current reactions, by Bliimlein and Tkabladze (3]. These calculations, based on the OPE, are 
extremely complicated, and we shall see presently how the same results can be obtained in a much simpler field- 
theoretic approach. 

The clue to this entire approach is contained in the classic paper of Ellis, Furmanski and Petronzio (EFP) which 
gave the first derivation of the dynamic higher twist corrections to unpolarized DIS in terms of amplitudes involving 
not just the "handbag" diagram of Fig.[Tl whose soft part is the quark-quark (qq) correlator $, but the higher order 
diagrams in Fig. [21 whose soft parts are the qqG and qqGG correlators respectively. 
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FIG. 1: The DIS "handbag" diagram involving the qq-correlator. 




FIG. 2: DIS Diagrams involving the qqG and ggGG-correlators. 

EFP begin with a brief discussion of a parton model, which they refer to as a "reference model", in which the 
active quark, momentum fc'^, emitted from the nucleon is on mass-shell. HandUng the kinematics exactly they arrive 
at expressions for the unpolarized structure functions i^i_2 in terms of the quark densities q{x) which are identical to 
those of Nachtmann . 

At first sight it seems surprising that such a naive model should give the exact results of the theory. But the point 
is - and this is something not stressed in the literature - that any off-shellness of the quark is a direct consequence 
of QCD, i.e. if the strong interaction coupling .g = then it follows that (taking the quark mass mq = 0) /c^ = [see 
Eq. ([5])]. Target-mass corrections are, by definition, kinematic in origin, therefore are independent of the value of g. 
Hence it is not miraculous that a model with = 0, i.e. equivalent to putting g — Q, should yield exact results for 
the target-mass corrections. However, it is crucial that the Lorentz structure is built into the model, as is done by 
EFP [8] . The implications of this are that the exact target-mass corrections must be derivable from the "handbag" 
diagram, Fig. [1] alone, since the diagrams of Fig. [5] vanish when (7 = 0. 

We shall carry out the derivation for polarized DIS and also show that when the Lorentz structure is respected the 
k"^ = {) "model" yields an unambiguous result for g2{x), namely the Wandzura-Wilczek (WW) result [l| 

5'2(a;Bj) = -.gi(2:Bj) + / dy , (1) 

where xbj is the well-known Bjorken variable for DIS. This suggests that the analytic continuation necessary in the 
OPE derivation of the WW result is in fact correct, and also implies that in a correctly formulated parton model, 
with on-shell quarks, g2{x) is exactly given by the WW expression. 

In this paper we shall also show that careful attention to the Lorentz structure in the parton model imposes strong 
constraints on the possible kr dependence of the so-called transverse momentum dependent (TMD) parton densities, 
that in recent years have received much emphasis (for an up-to-date review see e.g. Ref. Q and references therein). 
Indeed, taking into account the transverse momentum of partons is important for an understanding of the large 
transverse single spin asymmetries observed in many reactions. It is also essential in order to generate the parton 
orbital angular momentum which appears necessary as a consequence of the small contribution to the nucleon angular 
momentum provided by the parton spins. 

As we are going to show, imposing Lorentz covariance we are able to derive, for the unpolarized densities q{x, k"^), 
the longitudinal densities Ag(a;, fc|n) and the transversity densities /S.Tq{x, fey) or A'rpq{x, k'^) = hi{x, fc^), their 
dependence on kx from the functional form of the usual coUinear, fcr-integrated, parton densities q{x), Aq{x) and 
Axqix) [These relations are spelled out in detail in Eqs. ([57)l . (|6ip and (pD|) . and in Eqs. ([70)1 ]. Moreover, if 

perturbative evolution with is taken into account, then the evolution of the feT-dependent densities is entirely 
controlled via the known evolution of the coUinear densities. 

It should be noted that the functional form found in the parton model is quite different from the often used factorized 
form, typically f{x) e~^^T ^ though the latter may be a reasonable starting point for analyzing the presently available 
data. 
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Further, we obtain positivity and SofFer-type bounds [lO|, based on the fcy-dependent densities, which are stronger, 
in general, than those usuaUy imposed on the purely x-dependent, fc-r-integrated, coUinear densities, and we suggest 
that it would be interesting to impose these bounds on the collinear densities when extracting parton densities from 
deep inelastic and semi-inclusive deep inelastic scattering data. 

Therefore, as will become even more clear in the following, careful attention to the Lorentz structure in the parton 
model with on-shell partons, i.e. with = 0, has dramatic consequences. Given the standard belief that the soft 
functions appearing in QCD have support only in a very narrow range of values of around k^ — 0, it is tempting to 
suppose that the results derived in this paper may not be too different from those holding in a full QCD treatment. 

Covariant parton models have been discussed in the literature for many decades, beginning with the work of 
Landshoff, Polkinghorne and Short in 1971 [ll[, and Franklin in 1977 [l^l- These papers dealt only with the unpolarized 
structure functions ^"1^2(2;). The collinear spin-dependent structure functions 31.2(2;) were studied, in response to the 
"spin crisis in the parton model", by Jackson, Roberts and Ross |T^, who also commented upon some aspects of 
quark transverse momentum. More recenth^ a very detailed study of covariant models was initiated by Zavada in 
1997 [3], and later developed by Zavada [3, [lB| and Efremov, Teryaev and Zavada [l3|, to include the polarized 
parton densities as well as parton transverse momentum. While our paper was in its last phase of preparation we 
became aware of a further paper b y E fremov, Schweitzer, Teryaev and Zavada (which we shall refer to as ESTZ) 
and a very latest work by Zavada [l^ in which TMD parton distributions are discussed in detail. 

In these papers the analysis is based upon reasonable assumptions about the structure of the quark wave function 
in the nucleon rest frame. Because the treatment is covariant the results can then be translated to any other Lorentz 
frame. It is remarkable that many of the results of this kind of analysis are in complete agreement with our more 
general, frame independent, treatment. The principal differences are that: i) in our treatment the three leading 
twist densities, the unpolarized density q(x, the longitudinal density Aq{x, k'^) and the transversity density 
ATq{x, kx) or A'rpq{x, k'^) = hi{x, k'^) are independent, whereas in the Zavada et al. papers only two of these are 
independent; ii) we take into account target mass corrections. Because our approaches are so different, we believe it 
will be instructive to present our derivation in full. 

The plan of the paper is the following: in section [ll] we discuss the hadronic correlator for polarized DIS in the 
parton model, with g — 0, and give the general expression of the structure functions gi{x) and 52 (a;); in section Hill 
we evaluate explicitly gi (x) and 52 (x) including target-mass corrections and derive the WW relation; in section IIVI 
we extend our analysis to the TMD quark distributions and discuss the consequences of imposing Lorentz invariance; 
in section |V] we examine the positivity and Soffer bounds for the fcr-dependent unpolarized, longitudinally and 
transversely polarized quark distributions, discuss how their correct implementation leads to new, more stringent, 
bounds on the fcy-integrated collinear distributions and show some phenomenological implications of this approach; 
finally, in section IVTl we give some additional comments and conclusions. 



The correlator $y (P, 5;^), where P'^ and S*^ are the momentum and spin-polarization 4- vectors for the nucleon 
(5^ = —1) , fc^ is the quark 4-momentum and i,j are Dirac indices, is defined by 



To preserve colour gauge invariance a path-ordered Wilson line (or gauge-link) should be inserted between the quark 
fields. Since we are mostly interested in the g — case, we will take this operator as the identity. 

We ignore for the moment flavour and the quark charge - they are trivially reinstated at the end - and, as already 
stated, work with = 0. By partial integration we see that, when g = 0, 



II. THE HADRONIC CORRELATOR $ FOR g^O 




(2) 




(3) 



and that for any of the 16 Dirac matrices F (= /, 7'^, a^" , 7^75, 175) 



Tr[$r = . 



(4) 



Finally, again by partial integration, we get for all 




(5) 
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implying that 

$ cx (5(/c2) . (6) 

$ is a 4 X 4 matrix in Dirac spin space. The first attempt to write down its most general form was made by Ralston and 
Soper [2(j | who, however, missed one term j2lj |. The full dynamical twist- two structure, including the three so-called 
"naively T-odd" amplitudes, contains twelve terms [22,]. Eight of these, including all the "naively T-odd" amplitudes, 
are eliminated by the requirements of Eqs. ©-(O and two become related to each other. We are left with 

$ (P, S; fc) ^=%3 ^ + ^ (fc . 5) ^ ^5 + ^ a^- k, [{k ■ S) - {k ■ P) 5.] . (7) 

The scalar functions A^^^^n are, in general, functions of P • fc and k^. In our case, bearing in mind Eq. ©, we shall 
later use 

and 

9 ■ k 

^B.ii = - ^ (^) s{k') e[{p kf] , (9) 

where the factors are for later convenience and, following Rcf. 8], the 0- function ensures that the nucleon remnant 
has positive energy. Notice again the factor 6{k^) in these expressions. 

For the moment we shall only deal with the part relevant to gi^2{xBj) i-e. Ag. Polarized DIS is controlled by the 
antisymmetric part tV^" of the hadronic tensor, related to the structure functions .gi, 2(2^6]) via [23l | 



2M 



P-q 



W^r-^^'"'"^«^^/5ffi+ S^-^-^-^ gA, (10) 



P-q 



where q is the photon 4-momentum and our convention is €0123 — 1 



The contribution from the handbag diagram. Fig. [2 to W'f^ is then 



= e'^""" / d*fc(fcp + gp)(5[(fc + q)^]Tr(7,75<&) 



= — e 
4 

where 



^■^"^ / d^fc (fcp + qp) 6[{k + q)'] Tr(7, ^) 4? ' S) 



M 

e'-^-g, V , (11) 



V = j d^fc k, 5\{k + qf] As . (12) 

Note that the hadronic tensor W'^ in Eq. (fTTjl is electromagnetic gauge invariant. 

We introduce the standard two auxiliary null vectors and , = = 0, p • n = 1, in a slightly unusual way 
via 



1 M'^f 2F 

P'=D^^'--Q^'^'^ ^'^ = ;^('Z^ + eP''), (13) 



where ^ is the Nachtmann variable 



e= ; °^ (14) 

'l + 4M2a;2./Q2 



and 



i^^l + e ^^^^ 
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Note that 

q-p=^ q-n^-^ XBj = N=l-e. (16) 

Switching off the target-mass corrections thus corresponds to putting e = 0, D = N = 1. 

Since by definition, see Eq. (|12|) . I^a- is symmetric under /3 < — > a, its most general expression can be written as 

Ipa = M'^Bg gp„ + M'^Bpn ( ppn^ + p^np ) + Bp ppp„ + M^B^ npn^ . (17) 
On the other hand, from Eq. (fT2|) . bearing in mind the factor d{k^) in As, Eq. ([9|), and introducing the shorthand 

(X) =y d*kX6[{k + qf]As, (18) 

one finds that 

/n'^ Ip, = ( (A: • p) (fc • n) ) = M^Cpn 

P^'p'^Ipa = ((A:-p)') = M^Cp (19) 
n'^n" V = {(k-nf) = C'n . 

By comparing with Eq. (|17p and expressing the factors _Bi in terms of the Ct 's one can then write 

I/Sa = M'^Cpn [2{ppnc, + Pa-np) - gpcr] + C„ ppPa + M^Cp npn^ . (20) 
Moreover, from Eq. ^3]) . because of the factor e^^^^^qp S'^ in Eq. pTj) and the fact that P • = 0, we can replace 



in Eq. Thus 



iPa n„ ^ Pf3 ^ -^9/3 rif3 ^qp 



Ip, ^ - Af 2 Cp„ g^, + [ ANCpn -Cn+'ieCp]qpP, . (21) 

Putting this into Eq. (|11|) and comparing with Eq. (jlOp , we can read off directly expressions for gi + g2 and 52 and 
thereby obtain 

91 {xBi ) + 92 (xBj ) = 2 (P ■ g) Cp„ , (22) 

and 

51 (xBj ) = 2 (P • <z) { Cp„ - ^ [ 47V Cp„ - C„ + 4 e Cp ] } , (23) 

where we have used Eq. 

III. TARGET-MASS CORRECTIONS AND THE WANDZURA-WILCZEK RELATION 

In order to evaluate the coefhcients Cp„, Cp, C„ in Eqs. ((22|) . (p3)) we parameterize fc^ in the usual way, bearing in 
mind that — 0, 

kt' = xp'' + ^n>' + kl^, (24) 



where k'^ — —k^. Then 



d4fc = 1 / ^ d^fcj, dfc2 ^1 f^dk^ , (25) 
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and in terms of these variables 



k 



Cvn^{^,) = = ( ). (26) 



Bearing in mind Eqs. (fTB|) and one has 
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J[(fc + gf] = |<5[fc|.-Q2|(|-l)]. (27) 



We therefore integrate over fe|n, introduce 



via Eq. ([TT]) . and have 



[Dx-S,], (28) 



dx = ^l^ dv, (29) 

where the upper integration hmit comes from the 0- function in Eqs. ([8]), ([9]), and the lower limit from t] > x > In 
terms of rj the variables appearing in Eq. (|26p are given by 

We substitute the expressions for the coefficients, Eq. pS)) . into Eqs. P^ . (^5)) using Eq. ([HI), restore the quark charge, 
separate the leading twist terms, i.e. those that do not vanish when e ^ 0, and, inside the integral, order the rest 
according to the powers of e and rj involved, to obtain, for each flavour, 

ffi(^Bj) = dr, { e(-2e + ry) + e (8^' - ll^v + V) + ^'{-^e + n^V - ^v") + e'l-^?? + 2?7')} Ml) , (31) 

9iM+92{xBj)^el^ dv[^{^-7j) + eij{^-ij)}Mv)- (32) 
From here we also get a separate parton-model expression for g2'. 
92{xBi) = 1^ dv[a3^-2Tj) + e {^6e + mv - i^f) + e'(3e' - 10^'? + W) + £'(2^?? - 37?')} Mri) ■ (33) 

A. Neglecting target-mass corrections 

The expression for gi{xBj), Eq. (j3ip . must agree with the standard coUinear parton model result, with no target-mass 
corrections (NTM), in terms of the (longitudinal) polarized quark density, 

gf™M^'^Aq{xBi), (34) 
when target masses are neglected, i.e. when e — > and ^ — > ccbj. Thus 

Aq{xBj) = XBj / d?] {rj - 2 XBj) (fisiv) ■ (35) 



This relation between (fsifl) ^^d Aq{xBj) allows us to express the integrals in Eq. pip in terms of integrals over 
Aq{^), but also has powerful consequences for the structure of transverse momentum dependent parton densities, as 
will be discussed shortly. 



Turning to 5i(xBj) + <72(a;Bj), we have from Eq. ((32| : 

=2 /.I 



Bj 



Consider now [see Eq. ([55)) ] 



Changing the order of integration, we obtain < rj < rj' < 1) 



1 7 ^1 <.! 



Putting this relation in Eq. ([36)1 and using Eq. (l34l) yields 

NTM/ \ I NTM/ \ /""^ da;' NTM/ /\ 

5i (a;Bj)+52 (a^Bj)=/ — 5i (a;), 
which is the Wandzura-Wilczek relation in the absence of target-mass corrections. 

B. Inclusion of target-mass corrections 

Returning to Eqs. ^ we note first that from Eqs. §Ei), §B 



Then, using again Eq. ([35]) in the second step 



dr]r]ips{r]) = / dr] {r] - 2^ + 2^)ips{r]) 



Further, integrating by parts. 



d7] r]'^ ips{rj) ^5.1 drj 1] ipsi-q) + / dry / di]' rj' (ps{r]') 

^ dry ^ , , ^ dr] dry' 



Substituting Eqs. (|40| - ([42|) into Eqs. ([3T]l . (|32p we obtain expressions which, for reasons to be explained in 
subsection, we shall call the bare parton model (BPM) results 



„BPM 
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which, as can be shown after some algebra, agree exactly with the dynamical twist-two (i.e. taking c?„ = 0) OPE 
results of Ref . [6|] . Furthermore we find that 

C d r BPM/ N I BPM/ M ^ BPM/ n / Ar.\ 

? ^ [5i (a^Bj) + 92 (a^BjjJ = .9i (a^Bj) , (46) 
so that, bearing in mind Eqs. (|15|) and ([T6|) . 

<7f™(xBj) +5?™(.Bj) = /' ^5f™(^') , (47) 
and the Wandzura-Wilczek relation holds also when target-mass corrections are included Q. 



C. Target mass corrections and threshold behaviour 

There is a long-standing problem about the behaviour of the structure functions, both unpolarized and polarized, 
at XBj = 1. To see this most simply consider the first term in the expansion of Eq. in powers of A'P /Q'^, where 
we show the kinematic dependence in explicit detail 

5i(a^Bj)-^A<z[a; = e(xBj,Q2)]. (48) 

It is usually stated that one must have 5i(xBj = 1) = 0, implying the bizarre result that l^q{x) vanishes at 
X — ^(IjQ^) i.e. at a continuous infinity of dependent points. Of course without target mass corrections the 
vanishing of giix^i = 1) simply implies that lS.q(x = 1) = . 

There are no satisfactory prescriptions for avoiding this issue in the literature. Georgi and Politzer Q and Piccione 
and Ridolfi argue that higher twist terms must be taken into account in the region of large xbj , whereas Accardi 
and Melnitchouk [24i] impose some constraints on the virtuality of the struck quark. 

We would like to propose a very simple resolution to this problem. In the standard parton model treatment there 
is the underlying assumption that the struck quark and target fragments materialize into physical hadrons with 
probability one. This is tantamount to the assumption of equivalent completeness of partonic and hadronic states i.e. 

\parton) {parton\ = \hadron){hadron\ . (49) 

all parton all hadron 

states states 

While this might be reasonable for massless hadrons, it cannot be true when hadrons possess their physical masses. 
Thus what is needed is a non-perturbative function to express the failure of equivalent completeness. It has to 
express the fact that there is zero probability for the final partonic state (the struck parton and the target fragments) 
to hadronize if its total energy is too small to produce an inelastic hadronic event, and that the probability for 
hadronization is one when the energy is large enough. The simplest possibility is 9{xTh — a^Bj), where XTh is the 
maximum kinematically allowed value of xbj for inelastic scattering. We note that, strictly speaking, and excluding 
photon production which is higher order in the fine structure constant, this is not xxh = 1, but 

where fj, is the pion mass. 

Thus we propose that in a physical parton model (PPM) the BPM target mass dependent expressions, Eqs. P5)l . (|44p . 
for gf™ and gf^^ + gf^^ should be modified to 

g^^ixBi) = 5f™(^Bj) 0(xTh - XBj) (51) 



(52) 
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IV. RESULTS FOR THE TRANSVERSE MOMENTUM DEPENDENT QUARK DENSITIES 



In the standard coUinear parton model the transverse momentum kx is integrated over up to the large scale Q. 
But in the last few years many polarized reactions have been studied in a more general framework in which the 
fundamental parton densities q{x), Aq{x) and ATq{x) are allowed to depend on the intrinsic fc-r (see e.g. Refs. [ol. l25j 
and references therein). In all these approaches it is often assumed that the x and kr dependence can be factorized. 
This, as we shall see, is in contradiction with the Lorentz structure of the amplitudes, at least in the .g = case. We 
shall comment later on the more general situation. 

For unpolarized DIS the analogue of Eq. ([55)1 is P 



(53) 



q{xBj) = XBj / dri ifi{r]) . 
Moreover, on changing integration variables from rj to k^ via Eq. we see that q{x,k^) for which 

j d^kTq{x,k'^) = q{x) 

must be given by 



^T) = ^,V.{x + ^)e[x{l~x)M^~kl] 



xM^' 



(54) 



(55) 



Thus the x and k^ dependence are intimately linked. Moreover, the maximum size of Ifcrl is bounded and x-dependent, 
and in general fcf, < ~ 0.25 GeV^. 

Most importantly, if the functional form of q{x-Q^) is known (for consistency it must, as usual, satisfy q(l) ~ 0) then 
from Eq. ^ 



(56) 



Therefore we have the remarkable result that, in the g — case, the k^ dependence of q{x, fc|,) is completely 
determined by the xbj dependence of the standard coUinear quark density q{xBj)- 



q{x, k 



1 



ttA-P dx 



q(x) 



where, we remind the reader. 



■q — X + 



x—r] 



xM^ 



'[a;(l - x) Af 2 



krp] , 



(57) 



(58) 



Notice that Eqs. dMll-dMl) are in agreement with Eqs. (47), (48) of Ref. 

As a consequence, the average transverse momentum squared, ( fc^(x) )q, can be easily obtained: 



APx^ 
qix) 



(k'^{x))q — 



' q{y) 

dy 

y 



(59) 



For the case of polarized DIS the result is slightly more complicated. Given the functional form of Aq{xBj) we can 
obtain (psiv) '^ia Eq. (|40p. in agreement with Eq. (42) of Ref. ^16] : 



4{3Ag(r/) 



7? — Aq{Tj) + 2 
drj 



1 drj' 

n 



(60) 



The transverse momentum dependent polarized parton density is then given, via Eq. (j35p and the relation between 77 
and fc^ in Eq. jH]), by 



Aq{x, k'^) 



1 



7rAf2 



xM^ 



krp] , 



(61) 
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with 



J d^kr Aq{x, k^) = Aq{x) . (62) 

Just as for the unpolarized density, we see, via Eq. ([60| . that the dependence of Aq{x, k^) is completely determined 
by the XBj dependence of the coUinear polarized quark density Ag(a;Bj)- 

Of great interest at present is the transversely polarized quark density [l^] (also referred to as the transversity, see 
e.g. Ref. [IHl for a review) which is concerned with quarks transversely polarized along the y-direction inside a nucleon 
transversely polarized along OY: 

Arqix, kr) = gs^/Sy(x, kr) -q-sy/SviX' kr) 



V , k. 



Here we specify the parton distributions using the physically motivated probabilistic notation of Ref. (25j , referring 
to the polarization of a parton, whose spin direction (in its helicity frame whose axes are labelled ox, oy, oz (27j ) is 
indicated hy s = s^, Sy for transverse polarization, inside a nucleon moving along OZ in some fixed frame with axes 
labelled OX, OY, OZ (it could be the laboratory frame or the cm. frame of the 7*-nucleon collision), whose spin 
direction is indicated by S* = Sx, Sy for transverse polarization, namely 

Aqs/s{x, kr) = qs/six, kr) - q-s/s{x, kr) ■ (64) 

The connection with the more formal notation of Ref. [2^ will be given in section IIV Al (see also appendix C of 
Ref. [II])- In Eq- (ESI) 4> is the azimuthal angle of fc^ in the fixed frame. 
Note that 



/ 



d^feT At(7(x, fcr) = . (65) 



This is because Axqix, fcy) corresponds to having the quark polarized along oy perpendicular to its momentum and 
not along OY, perpendicular to the nuclcon's momentum. The A;T-dependent density which does integrate to the 
usual coUinear transversely polarized or transversity density is one in which the quark is polarized along the direction 
OY in the fixed frame: 

Ags^/Sy (a;, fc^) — A'rpq(x,k'^) in the notation of Ref. [IHl 

— hi{x,k^) in the notation of Ref. [12] (66) 

so that 



One finds 



Arqix) ^ <l^kThi{x,k^T)- (67) 



A'^qix, kl) = h,ix, kl) = (x + e[x{l - x) - k^] . (68) 



It follows that 



Arqix, kr) = — hi{x,k'^) cos (j) . (69) 



X 



Moreover, from Eq. (|67p one can show that 



1{Ha„,„-Aa„(,„}.-A[^]. ,™, 

Thus, via Eqs. (|68p . (|69p . the kr dependence of both hi{x, k'^) and Arqix, kr) is completely determined by the 
functional form of the standard coUinear transversity density Arq{x). 

It is an intriguing question whether, in the real case oi g ^ 0, it is possible to have a factorized form F{x)G{k'^) 
for the parton densities. Aside from the more complex tensorial structure when g ^ 0, the functional dependence will 
be controlled by expressions of the form 

dfc2$f k\x+^^^^^^ 



xM2 

where the range of fc^ is limited to a small region around k"^ = 0. While it seems unlikely that this could lead to a 
factorized form, we have been unable to demonstrate this mathematically. 
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A. Other transverse momentum dependent quark densities 



In this subsection we shall give expressions for several other TMD quark densities, which, while playing no role in 
the fcj^-integrated hadronic correlator for DIS, are important in oth er p rocesses, like semi-inclusive DIS or Drell-Yan 
dilepton production, and have been widely used in the literature [9l[20|. 

Related to Ag, we consider quarks polarized longitudinally (i.e. with helicity ±1/2) inside a transversely polarized 
nucleon (here it is irrelevant which transverse direction is involved). Then (sec also Ref. |22|) 



^(1sJSt{x, kr) =q+/ST{^: fcr) - q^/STi^, ^t) = — — — giT 



M 



with 



giT{x, fey) 



k 



For the quark polarized in the x-direction, we have 



(x, k^) I sin (j) ; 



ttIVP 



One also finds that 



2 



^ xlVP' 



For the function hiT{x,k'^) defined via 



hi{x, fc^) 



fe2 

hiT{x,kl) + 



h^j'{x,kT) 



we find 



hiT{x,k'^) 



ttAP 



xM^ 



ipii{x + e[x{l - x) M'-kl]. 



Finally, for a quark polarized along the a;-direction inside a longitudinally polarized nucleon. 



^<ls:,/Szi^^ f^r) = Qs^/+ix, k'^) - q_^j+{x, k^) 



\kT 



hii^{x, fey) , 



with 



Note that for consistency in the parton model, with g = 0, we have, in agreement with the results of EZTS, 

-r.2 



V 



hixix, k 



X 



hij^{x, kj,) 



(71) 



(72) 



(73) 
(74) 

(75) 
(76) 
(77) 
(78) 
(79) 



V. BOUNDS ON THE COLLINEAR DENSITIES 



In attempting to extract information on the coUinear polarized density Aq{x) and the transversity density Axq{x) 
by fitting data using parameterized forms for the densities, it is usual to impose the positivity bound 



on Aq{x), and the Soffer bound 



\Aq{x)\ <q{x) 



\ATq{x)\ < -[q{x) + Aq{x)] 



(80) 
(81) 



on Axq{x). However, it turns out that while these bounds are necessary they are not sufficient. 
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A. The new positivity bound 

The reason why Eq. (|80p is not sufficient is that from a probabihstic point of view one should expect 

|Aq(x,4)| <g(a;,4), (82) 

and while Eq. ((82|) implies Eq. ((80|l . the reverse is not true. Actually, as we shall now show, it turns out that it is 
sufficient to require just 

|A(7(a;,0)| <g(x,0). (83) 



Ag(a^, 0) = M^)0[x{l ~ x) Af2] , (84) 



g(x,O) = ^^3(x)0[x(l-x)M2]. (85) 

Since Eq. ([83)1 must hold for all x we find that (of course V3('7) > 0) 

< , 86) 

and this is sufficient to guarantee that Eq. (|82p holds. Thus, we suggest that in parameterizing the coUinear polarized 
density Aq{x) one should perhaps impose the stronger bound 



For then from Eq. ([61 



and from Eq. ([55] 



3 Aq{x) - X Aq{x) + 2 f ^ Aq{x') 



< q{x) ~ x qix) , (87) 
da; 



which follows from Eqs. ([86]) . ([60]) and ([56)1 and is in agreement with Eq. (48) of Ref. [TB|; of course, one must have 
q{x) — X q'{x) > . 

Notice that there is a further bound, originally derived in Ref. [lB|, Eq. (60), 

dx' 

Aq{x) + 2 — Aq{x') < q{x) , 
J X X 

which can be easily obtained from our Eqs. ([SO)) . ([34l) . (|53p and ([86| . 



B. Reinterpretation of the SofFer bound 

Soffer derived his bound on the collinear transversity density [l^ 

\ATqix)\<^[qix) + Aqix)], (89) 

by noticing the analogy between the quark correlator diagram and the diagram for the u-channel absorptive part of 
forward elastic quark-nucleon scattering. But collinear parton densities include all partons with transverse momentum 
fey — ^^'^ ^'^^s '^o^ consist of strictly forward moving partons. Hence, in the light of the existence of intrinsic 
kx, Soffer's analogy could be said to have been misinterpreted, and the correct bound should read 

\ATqix,0)\ = |A!r<z(x,0)| < ^ [g(x, 0) + Ag(a;, 0)] . (90) 

Now, from Eqs. ([69]) . ([68]) we have 

Arqix, 0) = A'^q{x, 0) = ^ ipnix) 0[xil - x) M^] , (91) 
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and substituting this and Eqs. (IM)) . (|85|) into Eq. ([90]) we obtain 

\<Pii{v)\ < ^ifsiv) ~ Vfsiv)] ■ (92) 

We suggest that in parameterizing models for the coUinear transversity density AtQ^x), the scalar functions 953,8,11 
should be calculated in terms of q{x), Aq{x) and Axqix), via Eqs. (|56|) . (|60)l . ((70)) . and, perhaps, should be forced to 
satisfy the bound in Eq. ([92]) . 

Finally, we have also explicitly verified that the bounds in Eqs. ([SS]). are equivalent to the bounds on the 
fcr-dependent functions derived in Ref. ^28(|, in the g = limit. 



C. Collinear and TMD quark distributions 



As a simple and illustrative example of our results, we consider here some of the consequences of imposing Lorentz 
invariance on the fex-dependence of TMD parton distributions. We will start from a well-known set of unpolarized 
and longitudinally polarized distributions extracted within the standard collinear approach, namely the GRV98 
and GRSV2000 ^30] sets respectively, and show the crucial role played by the new bounds. Since we are working 
in the g = approximation, we will consistently consider only u and d quark distributions, neglecting sea quarks. 
Moreover, since we are here mainly interested in the intrinsic kr dependence, we choose as reference scale a relatively 
low one, — 2 GeV^, in order to exclude large effects due to perturbative evolution. From these parameterizations 
for q{x,Q'^) and Aq{x,Q'^) we therefore generate, via Eq. ([57)1 and Eqs. ([SO)) . ([CT)) the fcr-dependent distributions 
q{x,k^,Q'^) and Aq{x,k'^,Q'^). 

Concerning transversity, we have to recall that at present no parameterization in the standard collinear approach is 
available. Double transverse spin asymmetries, Att, in the Drell-Yan processes p^^p^ £^£^+X, orp^^T ^ £+£^+X, 
are considered the best tool to this end, but are unfortunately still out of reach for present experimental setups. In 
the meantime, the study of transverse single spin asymmetries in polarized SIDIS processes in the TMD approach 
has been found to be the most promising way to get information on the unintegrated, fcT-dependent transversity 
distribution [3ll | . We will therefore utilize in our discussion the first ever parameterization available for the transversity 
of u and d quarks obtained by fitting, in the TMD approach, the so-called Collins azimuthal asymmetries measured 
in SIDIS pion and kaon production by the HERMES and COMPASS experiments. We have to stress, however, that 
this parameterization has been obtained by assuming a factorized x and fcT-dependence for the TMD distributions, 
which is at variance with our findings in the simplified, g = 0, limit. Given the still poor knowledge of the transversity 
distribution, we believe however that, despite this inconsistency, it is of interest to investigate the consequences of 




FIG. 3: Left panel: unpolarized (solid line), longitudinally polarized (dashed line) distributions and their ratio (dotted hne) 
for up quarks as a function of x at Ikr] = and — 2 GeV'^, obtained via Eq. (|57p and Eqs. ((60}, (|6ip starting from 
u(x) [GRV98] and Au{x) [GRSV2000]. Note the violation of the positivity bound, Eq. (|83]). Right panel: Soffer-type bound, 
11+ = (u + Au)/2, (solid line), transversity distribution (dashed line) and their ratio (dotted line) for up quarks as a function 
of a: at Ifcrl = and = 2 GeV^ obtained via Eqs. dTOjl, dH} starting from Atu{x) of Ref. [H. 



14 




FIG. 4: Same as in Fig.|3] but for down quarks. Note, in the right panel, the strong violation of the Soffer-type bound, Eq. (|90[). 



our approach. Starting from the most updated parameterization of the transversity distribution for valence quarks, 
Ref. [33|, and using Eqs. ^7U\i . (pB]) . we therefore generate the g = 0, fcr-dependent distribution A'j,q{x,k'^,Q'^). 

Our results are summarized in Figs. [3][71 In Fig. [31 left panel, we show the unpolarized and longitudinally polarized 
u quark distributions, and their ratio, as a function of x at fixed Ifc^l =0. Although the coUinear parton distributions 
adopted here fulfill the usual positivity bound, Eq. ((80|) . we see that the stronger positivity bound of Eq. ([83]) is indeed 
violated for some x values. This is most clearly seen by looking at the ratio Au/u (dotted line), where the violation 
reachs about 15%. In Fig. [3l right panel, we also show, again for u quarks, the transversity distribution Atu and the 
positive-helicity distribution u+ = (u + Au)/2 and their ratio (related to the Soffer bound), as a function of x at fixed 
Ifcrl = 0. In this case also the stronger, Soffer-type, bound, Eq. is fulfilled. 

Similar results are shown in Fig. 0] for d quarks. This time, while the stronger positivity bound derived in this 
section is always fulfilled, the Soffer-type bound, Eq. ([90)) . is instead violated over a large range of x values. 

In Fig. [n]we show, again for Au{x, k^) in the left panel and for A'rpu{x, fc|,) in the right panel, the \kT\ dependence 




FIG. 5: Left panel: unpolarized (solid line) and longitudinally polarized (dashed line) distributions for up quarks as a function 
of Iferl at X = 0.2 and = 2 GeV^, obtained via Eq. ^ and Eqs. ((6TJ starting from u{x) [GRV98] and Au{x) 

[GRSV2000]. Right panel: Soffer-type bound (solid line) and transversity distribution, Ay it, (dashed line) for up quarks, as a 
function of \kT\ at a:: = 0.2 and = 2 GeV^, obtained via Eqs. ((68]), starting from Atu{x) of Ref. [s^. 
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FIG. 6: Same as in Fig. [S] but for down quarks. 



of the u-quark distributions at a fixed value of x in the valence region, x = 0.2. For the longitudinal distribution, we 
see that the violation of the stronger positivity bound at |fc7-| = persists up to some relatively large values of Iferl- 
A similar behaviour should result for any value of x in the region of violation of the positivity bound at Ife^l = 0, 
shown in the left panel of Fig. [31 Notice also the node in Am, which is a general feature of Aq coming from Eq. (|6ip . 
Finally, notice from the right panel that, as expected, the validity of the Soffer-type bound at Ifcrl = guarantees its 
validity at any allowed value of Ifcrl- Analogous results for d quarks are presented in Fig. [6l and similar comments 
hold on the violation, this time, of the SofFer-type bound. 

For completeness, in Fig. [7] we show the resulting mean transverse momentum squared, {k'^{x))q, Eq. (|59p . as a 
function of x, for u and d quarks, at the same reference scale, = 2 GeV^. We briefly note that the shape of 
{k'^{x))q is quite reasonable and approximately flavour-independent, although its size is considerably smaller than 
expected from available phcnomenological analyses within the TMD approach [H, IsJ] . 

This is perhaps a hint that the larger values of ( fcy ) needed in the phcnomenological analyses may be due to hard 



0.06 




FIG. 7: Mean transverse momentum squared, {k^{x))q, Eq. (|59|l . as a function of x, for up (solid line) and down (dashed 
line) quarks at = 2 GeV^. 
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gluon emission and may therefore not be a measure of the size of the genuinely intrinsic transverse momentum. It 
would be interesting to construct models based upon such a picture. 



VI. CONCLUSIONS 



We have shown that the simple parton model has a very rich structure when due care is given to its Lorentz 
properties. On the basis that kinematic relations should not depend on the value of the strong coupling g, and thus 
should be valid when g = 0, we have been able to give a simple derivation of the exact target mass corrections for 
polarized DIS, without recourse to the operator product expansion, and have verified that the Wandzura-Wilczek 
relation holds when target mass corrections are included. We have also demonstrated that that there is a surprising 
and intimate connection between the intrinsic kx dependence of parton densities and the functional form of the fcr- 
integrated coUinear densities. Indeed the kx dependence can actually be derived from the coUinear densities. Using 
some published versions of the coUinear densities we have studied the implications of this relationship and found that 
our fcT^-dependent distributions, thus derived, are much narrower than those used in phenomenological work i.e. our 
( fcy } is much smaller. This suggests that some of the kx dependence usually described as intrinsic, may, in fact, be 
due to gluon radiation. 

One outcome of our work is the clear demonstration that the factorizcd form g{x)f{k'^) often used in phenomeno- 
logical studies, is untenable, at least when g = 0. While not proved, it seems likely that this result is quite general, 
though, as commented on, a factorized form does seem a reasonable starting point for the analysis of the presently 
available data. In addition, because of this interrelationship between kx and x dependence, we have been able to 
derive new positivity and Soffer-type bounds on the coUinear densities, which are stronger than the those usually 
utilized in DIS and SIDIS analyses. We have shown numerically that some of the published coUinear parton densi- 
ties, which satisfy the traditional bounds, do not always satisfy these stronger bounds. The particular case of the 
down-quark transversity may be of some interest. 
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